Divergent Photon Absorption in Topological Insulator Ultra-thin Films 
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We perform linear and non-linear photon absorption calculations in topological insulator ultra- 
thin films on a substrate. Due to the unique band structure of the coupled topological surface states, 
novel features are observed for suitable photon frequencies, including a divergent edge singularity in 
one-photon absorption process and a significantly enhancement in two-photon absorption process. 
The resonanct frequencies can be controlled by tuning the energy difference and coupling of the top 
and bottom surface states. Such unique linear and nonlinear optical properties make ultra-thin films 
of topological insulators promising material building blocks for tunable high-efficiency nanophotonic 
devices. 
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Time-reversal invariant topological insulators are new 
states of quantum matter with an insulating bulk state 
and gapless Dirac-type surface states [TH3]. A range 
of compounds have been found to be three-dimensional 
topological insulators [5HTU]. Layered Bi2Se3 is demon- 
strated to be a prototype three-dimensional topological 
insulator with a large insulating bulk gap of ~ 0.3 eV 
and metallic surface states with a single Dirac cone [5]- 
[7]- A thin layer of topological insulator is expected 
to be a promising material for high-performance opto- 
electronic devices such as photodetectors [11] and trans- 
parent electrodes [12] due to its spin-momentum-locked 
massless Dirac surface state, which is topologically pro- 
tected against time-reversal-invariant perturbations. 

Two-photon absorption (TPA) is a primary process of 
interest in various emergent photonics applications [T3l - 
[17] • For application purposes, a good TPA material must 
display large absorptive nonlinearities tuned within spe- 
cific spectral regions [18] [19]. To gain insight into the 
origin of large (degenerate) TPA coefficients fa we con- 
sider the expression for j3 in second order perturbation 
theory, which is proportional to the transition dipole mo- 
ments and joint density of states (JDOS): 
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where ip c , ipi and ip v are Bloch wavefunctions of the 
electrons in conduction, intermediate and valence bands, 
whose energies are £^ c (k c ), Ei(ki) and £? v (k v ). Hi is the 
interaction Hamiltonian and X is the light irradiance. In 
general, one can get large j3 when reaching the resonant 
condition {E\ — E v = huj). Moreover, sharp peaks in the 
frequency dependence of the TPA coefficient should occur 
at critical points of the JDOS, such as Van Hove singular- 
ities. Two-dimensional systems may offer a novel avenue 
for creating useful TPA materials as, unlike in three di- 



mensions, Van Hove singularities in two dimensions may 
induce divergent JDOS. 

In this work we show that topological insulator thin 
films could provide a powerful setting in which both lin- 
ear and nonlinear optical processes of interest are greatly 
enhanced and are also highly tunable. A key feature of 
topological insulators is the existence of robust topologi- 
cal surface states, in which electrons propagate as mass- 
less relativistic fermions [See Fig. [lji] . In an ultra-thin 
film (5nm or thinner for E^Sea) of topological insulator 
the top and bottom surface states are coupled, giving 
rise to an energy gap [7 , 20-22 . The coupling strength is 
controlled by the film thickness [20 and the energy differ- 
ence between the two surface states can be controlled by 
substrate or electrical gating. Such tunability makes the 
topological insulator thin film a unique two-dimensional 
electron system. Due to the coupling of surface states 
the conduction band minima and valence band maxima 
occur at the same (nonzero) wave- vector [Fig. ^p] , lead- 
ing to a divergent JDOS and thus a divergent one-photon 
absorption (OPA) at the gap frequency, illustrated by op- 
tical process a\ in Fig. ^p. Furthermore, by tuning the 
relative amplitude of the gap and the top-bottom energy 
difference, one can achieve a band structure in which a 
two-photon process (fa and fa in Fig. [TJd) is greatly en- 
hanced due to the existence of an intermediate band at 
the resonance frequency and the almost divergent JDOS 
of the initial and final states (5 in Fig. [TJd). With such 
properties, thin films of Bi2Se3 (and similar materials 
Bi2Te3 and Sb2Te3) are unique material building blocks 
for new nanophotonic devices. 

The low-energy physics of a topological insulator 
thin film is characterized by two copies of the topo- 
logical surface states on the top and bottom sur- 
faces. In the simplest topological insulators such as 
Bi 2 Se3, Bi 2 Te3, Sb2Te3, each surface state has a single 
Dirac cone. The surface state wavefunction is localized 
on the surface and decays exponentially away from the 
surface with a characteristic "penetration depth" £. For 
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FIG. 1. (color online) (a) Coupling of Dirac cones on oppo- 
site surfaces of a thin-film topological insulator; (b) Resulting 
surface band structure with divergences in joint density of 
states for l-(ai) and 2-photon (5) transitions. The possible 
1-photon optical transitions are indicated by arrows with an 
index of a±, ct2, «3, ct± and 2-photon by f3\, 02, 7, S. 

the Bi2Se3 family of materials £ ~ lnm. For ultra-thin 
films with thickness comparable with £, the overlap be- 
tween the surface-state wavefunctions from the two sur- 
faces of the film become non-negligible and hybridiza- 
tion between them has to be taken into account. In a 
thin film on a substrate, the chemical potentials of the 
top and bottom surfaces are inequivalent and the Dirac 
points are generically at different energies. Considering 
the inter-surface coupling and the chemical potential dif- 
ference one can write down the following low energy ef- 
fective model of the thin film which matches well with 
experiment [7( 123]: 

Ho = hvT z ®(cr x ky - a y k x ) + ^-T x ®l + A ih T z ®l, (2) 

where v is the Dirac velocity, and 5 and r are Pauli ma- 
trices acting on spin space and opposite surfaces, respec- 
tively. Time reversal invariance follows from [G, T~L] = 0, 
where = 10 i<j y K and K is complex conjugation. 
Ah is the hybridization between the two surface states. 
Aib is the inversion symmetry breaking, which can be 
substantially modified through electrical gating [7]. For 
simplicity we neglect the higher-order terms in k and we 
will discuss the effect of higher order terms at the end of 
the draft. The surface band dispersion is 

E(k) = Tv /(Hk|TA ib ) 2 + (A h /2) 2 , (3) 

The energy gap is E e( ± ge = Ah at the wavevector |k| = 
Afo/hv. Low-energy optical absorption by the surface 
states can occur with photon energy ranging from i? e dge 
to gap of the bulk bands. 

The electron-photon interaction is determined by min- 
imal coupling, i.e., by replacing k by k — eA/hc in the 
model Eq. ([2|, which leads to the interaction Hamilto- 
nian 

g 

Ui = —vr z <g> (cr x Ay - cr y A x ) , (4) 




FIG. 2. (color online) (a) Contour plot of OPA spectra for 
thin film topological insulators (logarithmic scale). Ah/ Aib 
can be tuned by film thickness and electrical gating. a(uj) is 
in units of i^al^ftZ. (b) Line cut for the 4QL Bi2Se3 thin 
film with Ah = 70 meV and Aib = 53 meV. The edge of the 
interband transition E e d ge is indicated by an arrow; important 
features are labeled ai —0*4. 

Here A = Ae is the optical vector potential with ampli- 
tude A and polarization e. The amplitude A is related 
to the light irradiance by X = s/e^uj 2 A 2 /2ttc. Here we 
neglect the small wave vector of light. Taking into ac- 
count the momentum conservation k for the initial \ip c ) 
and final \i/j v ) states, only vertical excitation processes 
contribute to absorption. The optical selection rules are 
ftHWil^i) oc k, (^ 4 |ft 3 |^M oc -k, ftMWil^i) ocinxk 
and (ip2 \T~Li 1^3 ) °t — ^ x k, where n is a unit vector 
normal to the momentum k. 

The energy spectrum of the effective model Q is 
shown in Fig. lb. For finite chemical potential offset Aib 
the coupling between the two surfaces leads to avoided 
crossing of the energy dispersion at finite wavevectors, 
and the valence band maxima and conduction band min- 
ima coincide on a ring |k| = A^/hv. This feature is 
essential for optical properties of the system, since the 
density of states of both conduction and valence bands 
diverge at the same wavevectors, enabling a divergence in 
the probability of the optical transition process marked 
by ol\ in Fig. lb. Other important optical transitions 
beside ol\ are also shown in Fig. lb. 

The OPA coefficient is 



where W\ is the transition probability for OPA per unit 

^i = |EE I W*l «i l^i)| 2 S (Efi(k) - M , (6) 
a f/i k 

and JS fi (k) = £f(k) - £|(k). Fig. 2a shows the OPA 
spectrum when the doping level is in the gap. The pro- 
cesses ai 5 2,3,4 contribute to different features at different 
frequencies as marked in Fig. 2a. In the following we will 
discuss the contribution of these processes in more detail. 

1. The onset of ql\ leads to the band edge singularity 
at the band gap energy £^ e dge- This singularity is 
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directly related to the JDOS divergence of the sur- 
face states in the form of p ex 1/ ^H 2 uj 2 — A|, due 
to the energy minimum at finite wavevector. For 
large photon frequencies, the E\ — >• E^ contribution 
to a(uj) is proportional to uj~ 3 . 

2. At the frequency Huj = yj A 2 b + Aj; , the transitions 
Ei — >> £4, £? 3 — >> ,£4 and £? 3 — » £? 2 occur at k = 0, 
as is labeled by 0L2 in Fig. lb. This process leads 
to a step discontinuity in the OPA spectra. OPA 
from E3 — >• E4 is exactly zero at y A 2 b + and 
has cj -2 dependence as Huj/A^ ^> 1. 

3. For frequency fkj ^> Aib, the transitions 0^3 (£?i — > 
,£4) and «4 (,£3 — >• £^2) occur far away from the 
avoided crossing wavevector A^/hv. In this limit 
the inter-surface coupling can be neglected, and the 
transition occurs within each surface. It has been 
studied in the graphene context (24J [25] that such 
a transition in a 2D Dirac fermion system leads 
to a universal frequency-independent contribution 

with a = e 2 /he the fine- structure con- 
stant. This contribution dominates the absorption 
probability in the high frequency limit as shown in 
Fig. 2b. 

For the surface state of a bulk topological insulator, 
the TPA coefficient is obtained by including all possible 
intermediate states in the surface bands, which leads to 
Ahick = (27T 2 /e u; co' 4 fi 3 ) {ye 2 1 c) . There is no resonant 
feature or Van Hove singularity. In a thin film new res- 
onant features appear due to the inter-surface coupling. 
There are two possible transitions E\ —> E2 — >• E4 (fa, 
fa) and E 3 ^ E 2 ^ E4 (7), as shown in Fig. ^p. Both 
transitions are included in the calculation of the TPA co- 
efficient if Huj > \/ A 2 h + (Ah/2) 2 , as both of them satisfy 
energy conservation. In Fig. [3] we show numerical results 
for TPA coefficients and corresponding optical processes. 
When A h /2 < Huj < ^A 2 h + (A h /2) 2 , the optical tran- 
sitions from valence bands to conduction band ,£4 is for- 
bidden by energy conservation, so that f3 = 0. As the 
photon frequency becomes larger, fauj) has two resonance 
frequencies, corresponding to transitions fa and fa, with 
the resonance condition E4 — E2 = E2 — Ei = Huj indi- 
cated by Eq. 0. These features represent large tunable 
absorptive nonlinearities, making thin film topological in- 
sulators promising TPA materials for applications. The 
double resonance is at = (5/i ± ^9A 2 h — 4A 2 l )/4/w, 
and it disappears when Ah /Aib reaches a critical value 
Ah/ Aib > 1-5 as shown in Fig. [3^l. In practice, the en- 
ergy Huj in Eq. needs to be replaced by Huj + iT in 
order to take into account the effect of carrier damping. 
Here, T is assumed to be a constant and inversely pro- 
portional to the dephasing time r. In our calculation 
we set r/A ib = 0.05. Fig. [3]d shows the TPA spec- 
trum for representative values of Ah/A^. fauj) shows 
as double resonance feature for Ah /Aib = 1.32, which 



corresponds to the experimental values observed for 4QL 
Bi 2 Se3 film [3 [23] . The strongest resonance feature oc- 
curs at fa since the optical transition matrix elements at 
fa are larger than at fa . In contrast there is no strong res- 
onance for Ah/Aib = 1.60. The process 7 does not satisfy 
the resonant condition for any photon frequency, so it has 
little contribution to the TPA coefficient. Fig. ^ shows 
the maximum of /3 (at fa) versus the parameter Ah/Aib- 
Although the resonant condition is satisfied at = Aib 
when Ah = 0, the associated transition from E\ — >> E2 
vanishes. The strongest j3 occurs at Ah/Aib ~ 0.5. 

The double resonance feature of the TPA coefficient 
from transition E\ — ^ E4 is due to the Rashba-type 
splitting, compared to the single resonance in bilayer 
graphene [26]. The Rashba splitting also gives rise to 
the divergent JDOS at the band edge, which has already 
been shown in the OPA coefficient. Obviously, there is 
no resonant intermediate states in the TPA process from 
Ei — > £?2, however, with the divergent JDOS and finite 
transition matrix elements at |k| = Aib, the TPA con- 
tributed by the process S around the gap is 

0(l->2)cx 1 (7) 

u 3 y (2hu) - Aj; 

It has a singular feature centered around Huj = Ah/2 due 
to the Van Hove edge singularity. It shows uj~ 3 depen- 
dence in the resonant region Huj ~ Ah/2 and uj~ 9 de- 
pendence in the off-resonant regions of Huj ^> A^, while 
TPA of gapless surface states in topological insulators 
has a uj~ 4 dependence for all photon frequency. 

Taking into account of both OPA and TPA, the change 
in the intensity of the light as it passes through the sam- 
ple is given by dX/dx = —aX—fiX 2 . The total absorption 
coefficient is given by a to tai = OL+f3Z. The nonlinearity is 
characterized by the ratio fal/a. which depends on the in- 
tensity X. Since for the TI film we have shown that /? has 
resonance features occurring at different frequencies from 
that of a, the ratio fa a can be greatly enhanced, enabling 
the realization of strong nonlinearity at low intensity of 
light. For the 4QL Bi 2 Se3, the parameters are estimated 
by r - 2.3 ps [27H22], T = 2.9 meV, A ib = 53 meV. 
For the resonance frequency of /3 1.25 < Huj/ A-^ < 1.8, 
the condition [3X/a ~ 1 can be satisfied for an (optical) 
electric field strength of 10 5 V/m. Such field strengths 
potentially could be realized at low incident optical pow- 
ers using photonic resonators with high ratio of quality- 
factor to mode- volume [30] . 

There are higher order terms such as hexagonal warp- 
ing term proportional to k 3 [31] in the surface state dis- 
persion relation. With such terms the gap due to avoid 
crossing of surface states is no longer uniform around 
the crossing wavevectors, and the JDOS becomes finite. 
However, the warping parameter is very small in Bi2Se3 
when the crossing energy is lower than 0.22 eV (defined 
respect to the Dirac point), so that the JDOS enhance- 
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FIG. 3. (color online) (a) Contour plot of TPA spectra of transitions from valence bands to E4 for thin film topological 
insulators (logarithmic scale). Ah/Aib can be tuned by film thickness and electrical gating; T/Aib = 0.05. f3 is in units of 
(47r 2 ^/e u; Af b )(ve 2 /c) 2 . (b) Line cut for A h /A ib = 1.32 (4QL Bi 2 Se 3 ) and A h /A ib = 1.6. A h /A ib = 1.32 has two resonances, 
labeled f3\ and #2- (c) The maximum of /3 vs. Ah/Aib . 



merit given by our low energy effective theory remains 
valid [32]. 

It worths mentioning that the gap of topological insu- 
lator thin films is always less than the bulk gap 0.3 eV, 
which is in the THz frequency range. However, Raman 
processes may lead to transition between E\ and £"2,4 in 
the optical frequency range where the intermediate states 
are high-energy bulk states [33, 34 , such Raman process 
should be greatly enhanced due to the divergent JDOS 
and could more conveniently be accessed in experiments 
and applications. 
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